Confinement effects on the vibrational properties of III-V and II- VI nanoclusters 
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We present a first-principles study of the confinement effects on the vibrational properties of 
thousand atoms (radii up to 16.2 A) colloidal III-V and II- VI nanoclusters. We describe how the 
molecular- type vibrations, such as surface-optical, surface-acoustic and coherent acoustic modes, 
coexist and interact with bulk-type vibrations, such as longitudinal and transverse acoustic and 
optical modes. We link vibrational properties to structural changes induced by the surface and 
highlight the qualitative difference between III-Vs and II- Vis. We describe the size dependence of 
the vibrations and find good agreement for Raman shifts and for the frequency of coherent acoustic 
modes with experiments. 
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Semiconductor nanoclustcr (NC), research is a rapidly 
growing field driven by the attractive idea to tailor ma- 
terial properties by acting on the morphology of the 
structures. The modification of the optical properties 
by merely changing the diameter of colloidal NCs (aka 
quantum dots) is one of the figurehead of nanostructure 
science. Indeed, the optical properties of NCs are well 
understood theoretically and well controlled experimen- 
tally. These experiments are usually performed at room 
temperature where vibrations are naturally involved. A 
solid understanding of the vibrational properties of NCs 
is a decisive step for the real world application of NCs, 
where the physical properties such as temperature broad- 
ening, loss of quantum coherence and relaxation of charge 
carriers are key components (6l-[l0j. However, despite a 
rather intense recent experimental work on the vibra- 
tional properties of NCs @, S [2^, 2^, ^ Hi] , the theo- 
retical understanding of the processes involved remains 
rather shallow. Most of the work has been carried out 
at the level of continuum models (e.g. Ref. [l|) or ef- 
fective potentials (e.g. Ref. [l-Q) that rely on parame- 
ters derived from the bulk. Only very recently, ab initio 
calculations of the vibrational properties have become 
possible for Silicon NCs with sizes comparable to the ex- 
perimental situation^. Unexpected effects of compet- 
ing bond-length contraction and undercoordination were 
unraveled 0. In general the study of this intermediate 
size regime, between bulk and molecule, is interesting 
not only from a technological point of view (how can 
we expect the properties to change?) but also from a 
fundamental point of view (how can molecular-type and 
bulk- type vibrations coexist?). 

In this work, we study the vibrational properties 
of InP, InAs, GaP, GaAs, CdS, CdSe, and CdTe 
NCs with radii ranging from 10.7 (Ini4iPi4oH|y2) to 
16.2 A (Cd36gS348H3QQ) via ab initio density functional 
theory (DFT). We find that structural relaxation effects 
are large and penetrate inside the core of the clusters for 
NCs up to radii of 15 A. These morphological changes are 
qualitatively different for Ill-Vs (where several surface 
layers contract inwards) and II- Vis (where long range 
ionic interactions lead to a large bond length distribu- 



tion at the surface). These structural effects dominate 
the vibrational properties and lead to a large blue shift 
with decreasing cluster size in III-Vs and strong broad- 
ening of the optical branches and a filling of the optic- 
acoustic phonon gap with surface states in II- Vis. We 
find a small red shift of Raman peaks for our largest 
cluster sizes, in agreement with experiment. We ana- 
lyze the bulk parentage of the cluster modes and find 
them to be fully mixed longitudinal-transverse and fully 
distinct acoustic-optical. We find two types of surface 
modes, the surface acoustic (SA) modes — as the lowest 
frequency modes — and the surface optical (SO) modes. 
We identify the coherent acoustic phonons and find good 
agreement with experiment and with results from the 
simple Lamb mo del. 

The NCs are constructed by cutting a sphere, centered 
on a cation, from the zinc blende bulk structure, and re- 
moving the surface atoms having only one nearest neigh- 
bor bond. The surface dangling bonds are terminated by 
pseudo-hydrogen atoms H* with a fractional charge of 
1/2, 3/4, 5/4, and 3/2 for group VI, V, HI and II atoms, 
respectively. The symmetry of the bulk is thereby 
conserved. 

The calculations are performed using the local density 
approximation (LDA), Trouiller-Martin norm-conserving 
pseudopotentials, and an ener gy c ut-off of 30 Ry for HI- V 
and 40 Ry for II- VI clusters [ll|. The forces are mini- 
mized to less than 3xl0~^ Ha/Bohr under constrained 
symmetry. The dynamical matrix elements are then cal- 
culated via finite difference and the vibrational eigen- 
modes and eigenvectors are obtained after diagonaliza- 
tion. More details can be found in [l^ . 

To analyze the eigenmodes in terms of bulk and surface 
contributions, we calculate the projection coefficients 
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where, Nc, Ng, and N are the core, surface, and total 
number of atoms, X^{I) represents the three components 
that belong to atom / from the 3Ai"-component eigenvec- 
tor. The sum of and is equal to one. We define 
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the surface atoms as the atoms belonging to the outer- 
most seven layers of the cluster (around 3 A thick). The 
vibrational eigenmodes arc further projected onto bulk 
eigenmodes M„_q as 
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Here, itj'j ^ are the 



with Un,q{I) = ul ^exp{iq ■ Ri) 

three components of the 6-component bulk eigenvectors 
corresponding to the atom of type t, phonon branch n 
and wave vector q. Rj is the atom position of atom / 
in the cluster. The bulk phonon modes used in Eq. ^ 
are calculated via density functional perturbation theory 

m- 

We have calculated the vibrational properties for a to- 
tal of 23 clusters made of InP, InAs, GaP, GaAs, 
CdS, CdSe, and CdTe. In this letter we extract the 
essence of these calculations and will select representative 
resuhs. The vibrational DOS of InP and CdS NCs and 
bulk are plotted with a broadening of 0.8 cm~^ in Fig.[T] 
Although formally transverse- and longitudinal acoustic 
(TA, LA), and optical (TO, LO) modes seize to exist in a 
cluster, the comparison with the bulk Fig. [IJe)(j) reveals 
an obvious bulk parentage. At this stage, we make three 
observations: First, the III-V materials [illustrated by 
InP in FigHIa-d)] show a blue shift of L0-, TO- and LA- 
derived cluster modes with decreasing size, while II- VI 
materials [illustrated by CdS in Figlljf-i)] show no such 
shifts. Second, the surface modes tend to completely fill 
the acoustic-optic phonon gap in II- Vis, while a gap re- 
mains for III- Vs. Third, the "broadening" of the bulk 
optical phonon branches induced by the confinement is 
larger for II- Vis than for III- Vs. 

In order to understand these three effects, we plot in 
Fig. m the nearest neighbor distances of relaxed IITVs 
and II- VI NCs as a function of the distance of the re- 
spective bond to the cluster center. In all cases, the 
bond length at the dot center is reduced through the 
presence of the surface. For all the NCs studied, be- 
sides CdSe and CdTe, wc find that for R=15 A the bond 
length at the center has almost recovered and is 99% of 
the bulk bond length. CdSe and especially CdTe recover 
their bulk bond length slower. The bond length distri- 
bution looks qualitatively different for the III-Vs and for 
the more ionic II- VI materials. For III-Vs the surface 
shells show a successive reduction of bond length, going 
outwards, while II- Vis show a large bond-length distri- 
bution, as previously observed for small clusters Q- The 
overall reduction of bond length in III-Vs along with the 
positive Griineisen parameters (describing the change in 
phonon frequencies with volume) explains the blue shift 
of the L0-, TO- and LA-derived cluster modes (first ob- 
servation above). The lack of shift in the TA modes 
stems from the small negative Griineisen parameter for 
this branch. We attribute the large bond-length distri- 
bution at the surface of IT Vis to the long-range ionic 
interaction and connect it to the broadening of optical 
branches and the filling of the phonon-gap (second and 
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FIG. 1: (Color online) Vibrational density of states (DOS) 
contributed by core atoms (black) and surface atoms (red) for 
(a)-(d) InP clusters, (e) bulk InP, (f)-(i) CdS clusters and (j) 
bulk CdS. 



third observations above). 

Experimentally, the vibrational properties are mostly 
explored via Raman scattering. To compare our results 
with experiment we have calculated the Raman intensi- 
ties in Fig. [3] using the theory proposed by Richter et al. 
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with our ab initio projection coefficient ^ and dot 
vibrational frequency w". For the natural Lorentzian 
linewidth we use 2 cm~^ and the coefficients C^ ^ are 
summed up to Aq w 1/ (2i?) . Although NCs with diam- 
eter down to 2 nm can be produced, Raman shifts have 
been reported for NCs with more than 3 nm diameter. 
Only our largest NCs reach this size and we focus in Fig. 
[3] on clusters with radii around 15 A and 16 A. For CdS 
and CdSe we obtain a small red shift of the most intense 
Raman peak with decreasing NC size. This is in agree- 
ment with experiments (CdS [1^1, CdSe [13): although 
experiments are performed on ensemble of NCs and a 
proper comparison would require an average over several 
NCs, allowing for small variations in shape and surface 
morphology. A nearly constant Raman peak [l^ [l^ was 
reported for large InP clusters, which is consistent with 
our results for radii from 14.5 to 15.6 A in Fig. [T](c) and 
(d). Wc note that 4 meV blue shifts have been reported 
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FIG. 2: (Color online) Bond length distribution as a function 
of their distance to the dot center. For (a) InP, (b) CdS, (c) 
III-Vs and (d) II- Vis. LDA bulk bond lengths are given as 
dashed lines. 
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FIG. 3: Calculated Raman spectra for CdS and CdSe 
nanocrystals for different sizes. 



for InP (110) surfaces [201 ■ 

We now proceed to quantitatively analyze the bulk 
parentage of the cluster modes according to the pro- 
jection in Eq. From Fig [T] we can see that besides 
some modes in the gap and at very low frequency, the 
modes are generally mixed bulk-surface (red and black) 
modes. This is true up to our largest clusters, although 
the surface character of "bulk" modes has significantly 
decreased compared to the small clusters. In Fig. 0] (b) 
and (d), we show the vibrational DOS of InP and CdS 
clusters, where the character (TA, LA, TO, and LO) of 
the modes are color coded (black, red, green, and blue, 
respectively). The phonon dispersion curves of bulk InP 
and CdS are plotted in Fig. g] (a) and (c). From Fig. H 
(b) and (d), we see that the cluster eigenmodes have ei- 
ther optical or acoustic character with virtually no mix- 
ing. Interestingly, the "pure" surface modes which vi- 
brate with frequencies below the lowest core modes have 
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FIG. 4: (Color online) Bulk phonon dispersion curves for (a) 
InP and (c) CdSe. Vibrational DOS projected onto the bulk 
modes for (b) InP and (d) CdS NCs. 



purely acoustic character and the modes within the gap, 
between the acoustic and optical modes, have purely op- 
tical character. These surface modes can be described 
as "surface acoustic" (SA) and "surface optical" (SO) 
modes, respectively. We note that in principle, a mixing 
of acoustic and optical modes is allowed from symmetry 
arguments. Indeed, both acoustic and optical modes have 
the same symmetry within the point group[2l|, and we 
can see small mixing around the frequency of SO modes 
in Fig.|4](b) and (d). From Fig.|4l we also notice that the 
transverse and longitudinal character of the mode gets 
lost in the cluster, which shows fully mixed transverse- 
longitudinal modes. This mixing is introduced by the 
lack of translational symmetry. In NCs, the propagat- 
ing waves with various wave vector q arc combined into 
standing waves, pictorially via reflections at the bound- 
ary. This result is consistent with the experimental mea- 
surements p^] . where Raman activity was measured for 
frequencies around the TO-branch. 

Besides the mixing of transverse and longitudinal char- 
acters, confinement has a direct effect on the lowest non- 
zero vibrational frequency. Similarly to classical waves 
confined in one dimension, the vibrational frequency of 
the longest wavelength in a NCs can be described as 
/miri = vj^d. Here, v is the velocity of sound and d 
is the diameter of the cluster. In Fig. [5l we plot the size- 
dependent lowest frequencies /mm calculated from the 
longitudinal and transverse sound velocities as solid and 
dashed curves. As circles and crosses wc plot the low- 
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FIG. 5: (Color online) Lowest modes with bulk character 
(circles), surface acoustic modes (crosses), lowest breathing 
modes (triangles) and experimental results [2^, [SJ] for the 
coherent acoustic modes (square and diamonds). Lowest 
spheroidal mode according to the Lamb model (dashed dotted 
line), according to the confined bulk model using the sound 
velocity of the TA- (solid line) and LA-branch (dashed lines). 



est core and SA modes obtained from the DFT calcula- 
tions. We see that the lowest modes with dominant bulk 
character follow closely the analytic 1/R dependence. 
The SA modes have the lowest frequency since they can 
run over the surface of the cluster. They are however 
strongly affected by the morphology of the surface and 
their frequencies are not monotonous with cluster size. 
One last important type of modes are the so called coher- 
ent acoustic modes. These are breathing modes, where 
all the atoms vibrate in phase. These modes have been 
observed with Raman spectroscopy, far-infrared absorp- 
tion spectroscopy, and resonant high-resolution photolu- 
minesccncc spectroscopy [2^-[2^, and are now the center 
of attention when the manipulation of spins and the spin- 
decoherence is investigated. We plot our results as trian- 
gles in FiglU the experimental results as diamonds and 
the results of the Lamb modeljl^l as dashed dotted line. 
Our results are in good agreement with the experimen- 
tal results (although our clusters sizes are still somewhat 
smaller than experiment in the case of InAs) and with 
the simple Lamb model. In the continuum Lamb model, 
only the bulk transverse and longitudinal sound velocity 
enters. It may therefore be surprising to give such a good 



description, although we described that the clusters un- 
dergo strong lattice relaxations, absent in a continuum 
description. We attribute the validity of Lamb model to 
the small dependence of the sound velocities on the bond- 
length (2% change for a lattice deformation of 0.1 A) and 
the collective nature of the modes, where an averaging of 
the lattice deformations (especially for II- Vis) is a good 
description. 

By studying the vibrational properties of III-V and 
II- VI NCs with up to one thousand atoms via ab initio 
DFT we identify the following confinement effects. The 
LA. TO and LO-derived cluster modes of III-V clusters 
significantly blue shift with decreasing cluster size. For 
II- VI clusters this shift is absent but the broadening of 
bulk derived modes is significant and the gap between 
optical and acoustic phonons is filled by surface modes. 
We can clearly ascribe these observations to the large re- 
laxation of the clusters dominated by: an inward relax- 
ation of the surface penetrating deep inside the cluster in 
case of the III-Vs and a large distribution of bond length 
at the surface of II- Vis. These strong confinement ef- 
fects tend to disappear for clusters with more than 1000 
atoms, where a small red shift of the Raman peaks (we 
can see for our largest clusters) remains, due to a soften- 
ing in response to undercoordination. The modes are all 
fully transverse-longitudinal mixed but surprisingly dis- 
tinct in terms of their optical acoustic characters. We 
find surface optical modes in the phonon gap and sur- 
face acoustic modes as the lowest frequency modes. The 
coherent acoustic phonons are identified and found to be 
in good agreement with results from the Lamb model 
and experiment. We explain why the simple model by 
Lamb gives an accurate description in case of the breath- 
ing modes while the vibrational properties of small NCs 
are poorly described by continuum models in general. 
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